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& ■ Abstract 

o : 

^-H , In this paper, we develop a general theory on the coverage probability of random inter- 

vals defined in terms of discrete random variables with continuous parameter spaces. The 
£-H ' theory shows that the minimum coverage probabilities of random intervals with respect to 

m • corresponding parameters are achieved at discrete finite sets and that the coverage probabili- 

. ties are continuous and unimodal when parameters are varying in between interval endpoints. 

\ The theory applies to common important discrete random variables including binomial vari- 

^ able, Poisson variable, negative binomial variable and hypergeometrical random variable. The 

theory can be used to make relevant statistical inference more rigorous and less conservative. 

o : 

T— I ■ 

^ ! 1 Binomial Random Intervals 

Let X be a Bernoulli random variable denned in a probability space (0, & ', Pr) such that Pr{X = 
; 1} = p and Pi{X = 0} = 1 — p where p € (0, 1). Let X±, ■ ■ ■ , X n be n identical and independent 

samples of X. In many applications, it is important to construct a confidence interval (L, U) such 
that Pr{L <p<U\p}^l-5 with 5 G (0, 1). Here L = L(n,S,K) and U = U(n,S,K) are 
multivariate functions of n, 5 and random variable K = Y^l=\ -^i- To simply notations, we drop 
the arguments and write L = L(K) and U = U (K). Also, we use notation PrjL(A') < p < U (K) \ 
p} to represent the probability when the binomial parameter assumes value p. Such notation is 
used in a similar way throughout this paper. We would thus advise the reader to distinguish this 
notation from conventional notation of conditional probability. 

Clearly, the construction of confidence interval is independent of the binomial parameter p. 
But, for fixed n and S, the quantity Pr{L(A') < p < U(K) \ p} is a function of p and is 
conventionally referred to as the coverage probability. In many situations, it is desirable to know 
what is the worst-case coverage probability for p belonging to interval [a, b] C [0,1]. For this 
purpose, we have 



X 



*The author had been previously working with Louisiana State University at Baton Rouge, LA 70803, USA, 
and is now with Department of Electrical Engineering, Southern University and A&M College, Baton Rouge, LA 
70813, USA; Email: chenxinjia@gmail.com 



1 



Theorem 1 Suppose that both L(k) and U(k) are monotone functions ofk £ {0, 1, ■ • ■ , n}. Then, 
the minimum ofPr{L(K) < p < U(K) \ p} with respect to p G [a, b] is attained at the discrete set 
{a, b} U {L(k) 6 (a, 6) : < k < re} U {U(k) G (a, 6) : <k < n}. 

We would like to emphasis that the only assumption in Theorem 1 is that both L(k) and 
U(k) are either non-decreasing or non-increasing with respect to k. The interval (L(K),U(K)) 
can be general random interval without being restricted to the context of confidence intervals. 
This theorem can be generalized as Theorem 7 in Section 4. The application of the theorem is 
discussed in the full version of our paper [3]. Specially, Theorem 1 can be applied to the sample 
size problems studied in [T]. 

For closed confidence interval [L, U], it is interesting to compute the infimum of Pi{L(K) < 
p < U (K) | p} with respect to p G [a, b] C [0, 1]. For this purpose, we have 

Theorem 2 Suppose that both L(k) and U(k) are monotone functions ofk G {0, 1, • • • , n}. Then, 
the infimum of Pr{L(K) < p < U (K) \ p} with respect to p G [a, b] equals the minimum of the set 
{C(a), C(b)} U {Cuip) : p G J3u} U {C L (p) : p G J2 L }, where 

B v = {U(k) e(a,b):0<k< n}, J3l = {L(k) £ (a, 6) : < k < n}, C(p) = Pt{L(K) <p< U{K) \ p}, 
Cuip) = Pr{L(K) <p< U{K) \ P } and C L {p) = Px{L{K) <p< U{K) \ p}. 

It should be noted that the only assumption in the above theorem is that both L(k) and U (k) 
are either non-decreasing or non-increasing with respect to k. The interval [L(K), U(K)] can 
be general random interval without being restricted to the context of confidence intervals. This 
theorem can be considered as a specialized result of Theorem 7 in Section 4. 

2 Poisson Random Intervals 

Let X be a Poisson random variable defined in a probability space (f2, J^, Pr) such that 

X k e~ x 

Pr{X = k} = — — — , k = 0,1,2,- ■■ 

where A > is called the Poisson parameter. Let X±,--- ,X n be n identical and independent 
samples of X. It is a frequent problem to construct a confidence interval (L, U) such that Pr{L < 
A<t r |A}«l — 5 with 5 G (0,1). Here L = L(n,5,K) and U = U(n,5,K) are multivariate 
functions of re, 5 and random variable K = X^ILi^i- For simplicity of notations, we drop the 
arguments and write L = L(K) and U = U(K). For fixed n and 5, the coverage probability 
Pr{L(K) < A < U (K) | A} is a function of A. The worst-case coverage probability with respect 
to A belonging to interval [a, b] C [0, oo) can be obtained by the following theorem. 

Theorem 3 Suppose that both L{k) and U{k) are monotone functions of non-negative integer k. 
Then, the minimum of Pr{L(i^) < A < U(K) \ A} with respect to A G [a, b] is attained at the 
discrete set {a, b} U {L(k) G (a, b) : k > 0} U {U(k) G (a, b) : k > 0}. 
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It should be emphasized that the interval (L(K), U(K)) can be general random interval with- 
out being restricted to the context of confidence intervals. The only assumption in the above 
theorem is that both L(k) and U(k) are either non-decreasing or non-increasing with respect to 
k. This theorem can be generalized as Theorem 7 in Section 4. The application of the theorem is 
discussed in the full version of our paper [3] for the sample size problems studied in [2] . 

For the exact computation of the infimum of coverage probability Pt{L(K) < A < U(K) | A} 
for the closed confidence interval [L, U], we have 

Theorem 4 Suppose that both L(k) and U(k) are monotone functions of non-negative integer k. 
Then, the infimum ofPv{L(K) < A < U(K) | A} with respect to A £ [a, b] equals the minimum of 
the set {C{a), C{b)} U {Q/(A) : A £ £ v } U {C L (A) : A £ £ L } where 

B v = {U(k) e(a,b):k> 0}, £ L = {L(k) e(a,b):k> 0}, C(A) = Pv{L(K) < A < U(K) | A}, 
Cu(X) = Pr{L{K) < A < U{K) | A}, C L (X) = Pi{L(K) < A < U(K) \ A}. 

In Theorem 4, the interval [L(K), U(K)} can be general random interval without being re- 
stricted to the context of confidence intervals. This theorem is a special case of Theorem 7 in 
Section 4. 

3 Negative-Binomial Random Intervals 

Let K be a negative binomial random variable such that 

Pr{K = k}=(^ + r k ~ *V (1 - P) k , k = 0, 1, • • ■ 

with parameter p £ (0, 1) and r > 0. In the special case that r = 1, a negative binomial random 
variable becomes a geometrical random variable. For the coverage probability of open random 
interval (L(K), U(K)) for a negative binomial random variable K, we have 

Theorem 5 Suppose that both L(k) and U(k) are monotone functions of non-negative integer k. 
Then, the minimum ofPr{L(K) < p < U(K) \ p} with respect to p £ [a, b] C (0, 1) is attained at 
the discrete set {a, b} U {L(k) £ (a, b) : k > 0} U {U(k) £ (a, b) : k > 0}. 

This theorem can be readily obtained by applying Theorem 7 of Section 4. For the coverage 
probability of closed random interval [L(K), U(K)] for a negative binomial random variable K, 
we have 

Theorem 6 Suppose that both L(k) and U(k) are monotone functions of non-negative integer 
k. Then, the infimum of Pr{L(K) < p < U(K) \ p} with respect to p £ [a, b] C (0,1) equals the 
minimum of the set {C(a), C(b)} U {Cu(p) : p £ £?u} U {Ci{p) : p £ J3l}, where 

B v = {U(k) e(a,b):k> 0}, i2 L = {L{k) e (a,b) : k > 0}, C(p) = Pi{L(K) <p< U(K) | p}, 
C v (p) = Pt{L(K) <p< U{K) I p} and C L (p) = Pv{L{K) <p< U(K) \ p}. 
This theorem can be easily deduced from Theorem 7 of next section. 
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4 Fundamental Theorem of Random Intervals 

In previous sections, we discuss coverage probability of random intervals for specific random 
variables. Actually, the results can be generalized to a large class of discrete random variables. In 
this direction, we have recently established in [3] the following fundamental theorem of random 
intervals. 

Theorem 7 Let K be a discrete integer-valued random variable parameterized by 8 G 0. Let 
[a,b] be a subset of 0. Let L(K) and U(K) be functions of random variable K. Let J2u denote 
the intersection of [a, b] and the support of U (K). Let £Il denote the intersection of [a, b] and the 
support of L(K). Suppose the following assumptions are satisfied. 

(i) For any integers k and I, Pr{/c < K < I | 9} is a continuous and unimodal function of 
9 G [a, b]. 

(ii) For any 9 G [a, b\, there exist intervals Lq and I' e of real numbers such that {L(K) < 9 < 
U(K)} = {K G I g } and {L(K) <9< U{K)} = {K G I' g }. 

Then, the following statements hold true. 

(I) The minimum of Pr{L(K) < 9 < U(K) \ 9} with respect to 9 G [a, b] is attained at the 
discrete set {a, b} U £?u U 

(II) The infimum ofVi{L(K) < 9 < U{K) \ 9} with respect to 9 G [a, b] equals the minimum 
of the set {C(a), C{b)} U {C v {9) : 9 G £u} U {C L {9) : 9 G £ L }, where C{9) = Yi{L{K) <9< 
U(K) | 9}, C v {9) = Yt{L{K) < 9 < U{K) \ 9} and C L {9) = Pr{L(K) <9< U(K) \ 9}. 

(III) For both open and closed random intervals, the coverage probability is continuous and 
unimodal for 9 G (9',9"), where 9' and 9" are any two consecutive distinct elements of {a, b} U 

This theorem is proved in Appendices A. The notion of unimodal functions used in Theorem 
[7] is described as follows: 

A function of 9 is said to be unimodal for 9 G [a, b] if there exists 9* G [a, b] such that the 
function is non-decreasing for 9 G [a, 9*] and non-increasing for 9 G [9*, b]. 

It should be noted that the assumptions (i) and (ii) are satisfied for the following common 
discrete random variables: 

• Binomial random variable; 

• Poisson random variable; 

• Geometrical random variable; 

• Negative binomial random variable. 

A sufficient but not necessary condition to guarantee assumption (ii) is that L(.) and U(.) are 
monotone functions. 
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5 Hypergeometrical Random Intervals 

So far what we have addressed are random intervals of variables with continuous parameter spaces. 
In this section, we shall consider random intervals when the parameter space is discrete. We focus 
on the important hypergeometrical random variable. 

Consider a finite population of N units, among which M units have a certain attribute. Let 
K be the number of units found to have the attribute in a sample of n units obtained by sampling 
without replacement. The number K is known to be a random variable of hypergeometrical 
distribution. 

It is a basic problem to construct a confidence interval (L, U) with L = L(N, n, 5, K) and 
U = U(N, n, 8, K) such that Pr{L < M < U \ M} ~ 1 — 5. Here, U and L only assume integer 
values. For notational simplicity, we write L = L(K) and U = U(K). In practice, it is useful 
to know the minimum of coverage probability Pr{L < M < U \ M} with respect to M € [a,b], 
where a and b are integers taken values in between and N. For this purpose, we have 

Theorem 8 Suppose that L(0) < L(l) < • • • < L(n) and U(0) < U(l) < < U(n). Then, the 
minimum of Pr{L(K) < M < U(K) \ M} with respect to M € [a, b] is attained at the discrete 
set I UL , where I UL = {a,b} U {L(k) € (a, b) : < k < n} U {U(k) € (a, b) : < k < n}. 
Moreover, Pi{L(K) < M < U(K) \ M} is unimodal with respect to M in between consecutive 
distinct elements of Ijjl- 

For a proof, see Appendix B. In Theorem[8l the interval (L(K), U(K)) can be general random 
interval without being restricted to the context of confidence intervals. This theorem can be 
applied to the sample size problems discussed in [3] . 

A Proof of Theorem 7 

We need some preliminary results. 

Lemma 1 Suppose that {6' < L(K) < 6"} = {9' < U(K) < 0"} = 0. Then, 

{L(K) <6< U{K)} = {L{K) <9< U(K)} = {L(K) < 6' < U(K)} = {L(K) < 6" < U(K)} 
for any 9 G (9', 6"). 

Proof. By the assumption of the lemma, we have {L(K) < 9"} = {L(K) < 9'} U {9' < L(K) < 
9"} = {L{K) < 9'} and {9' < L(K) < 9} C {9' < L(K) < 9} C {$' < L{K) < 9"} = for any 
9 € (9', 9"). Consequently, 

{L{K) <9} = {L{K) < 9'} U {9' < L(K) < 9} = {L{K) < 9'} = {L(K) < 9"}, (1) 
{L{K) <9} = {L(K) < 9'} U {9' < L(K) < 9} = {L(K) < 9'} (2) 
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for any 9 G (0', 9"). Combining (TJ]) and ^ yields 

{L{K) < 0} = {L(K) <9} = {L{K) < 9'} = {L(K) < 9"}, V0 G {9', 9"). (3) 

By the assumption of the lemma, we have {U{K) > 9'} = {U(K) > 9"} U {9' < U(K) < 
9"} = {U(K) > 9"} and {9 < U(K) < 9"} C {9 < U(K) < 9"} C {9' < U(K) < 9"} = for any 
9 G (0', 6"). Consequently, 

{U(K) > 0} = {U{K) > 9"} U {0 < U(K) < 9"} = {U{K) > 9"} = {U(K) > 9'}, (4) 
{U{K) >9} = {U(K) > 9"} U {0 < U(K) < 9"} = {U{K) > 9"} (5) 

for any 9 G (9', 9"). Combining fll]) and ([5]) yields 

{U{K) >9} = {U(K) >9} = {U{K) > 9'} = {U{K) > 9"}, W G {9', 9"). (6) 

By taking intersection of events and making use of ([3]) and Q, we have 

{L(K) <9< U{K)} = {L{K) <9< U(K)} = {L(K) < 9' < U(K)} = {L(K) < 9" < U(K)} 

for any 9 G (9' , 9"). This completes the proof of the lemma. □ 

Now we are in a position to prove Theorem 7. First, we shall show statement (I). Let 9' < 9" 
be two consecutive distinct elements of {a, b} U £} v U Then, {9' < L(K) < 9"} = {9' < 
U(K) < 9"} = and by Lemma HJ we have 

Pr{L(iT) < 9 < U{K) \ 9} = Pv{L(K) < 9' < U{K) \ 9} = Pv{L(K) < 9" < U{K) | 9} 

for any 9 G (9', 9"). By assumption (ii), for any 9 G [a, b], event {L(K) < 9 < U(K)} can be 
expressed as an event that K is included in an interval. This implies that {L(K) < 9" < U(K)} 
can be expressed as an event that K is included in an interval. It follows from assumptions (i) 
that Pr{L(J^) < 9" < U(K) | 9} is continuous and unimodal function of 9 G [a, b]. Hence, for 
9 G {9\ 9"), letting < rj < min(0 - 9', 9" - 9, ^1), we have 9' + rj < 9 < 9" - rj and 

Pi{L(K) < 9" < U(K) | 9} > mm (Pr{L(K) < 9" < U{K) \ 9' + rj}, Vv{L{K) < 9" < U(K) | 9" - r)}) . 
By virtue of the continuity of Pr{L(K) < 9" < U (K) \ 9} with respect to 9 G [a, b], we have 
limPr{L(i<Q < 9" < U(K) \0' + rj} = Pv{L(K) < 9" < U(K) | 9'}, 

TjlO 

limPr{L(lf) < 9" < U(K) | 9" -r}} = Pv{L(K) < 9" < U(K) | 9"}. 

Therefore, 

Pr{L(K) < 9 < U{K) \ 9} = Pv{L(K) < 9" < U(K) | 9} 

> min {Pt{L(K) < 9" < U(K) | 9'}, Pr{L{K) < 9" < U{K) \ 9"}) 
= min (Pt{L(K) < 9' < U{K) \ 9'}, Pt{L(K) < 9" < U{K) \ 9"}) 

> min (Pt{L(K) < 9' < U(K) \ 9'}, Pt{L(K) < 9" < U(K) | 9"}) 
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for any 9 G (9', 9"). This implies that the minimum of Pi{L(K) < 9 < U(K) \ 9} with respect to 
9 G [9', 9"] is achieved at either 9' or 9". It follows that statement (I) is true. 

Next, we shall show statement (II). Let 9' < 9" be two consecutive distinct elements of 
{a,b} U J% U B L . Then, {& < L(K) < 9"} = {9' < U{K) < 9"} = and, by Lemma QJ we 
have Px{L(K) < 9 < U(K) | 9} = Pr{L(K) < 9' < U{K) \ 9} = Pr{L(K) < 9" < U(K) \ 9} 
for any 9 G (9', 9"). By assumption (ii), for any 9 G [a, b], event {L{K) < 9 < U(K)} can be 
expressed as an event that K is included in an interval. This implies that {L{K) < 9" < U(K)} 
can be expressed as an event that K is included in an interval. It follows from assumptions (i) 
that Pr{L{K) < 9" < U(K) \ 9} is continuous and unimodal function of 9 G [a,b]. In the course 
of proving statement (I), we have shown that 

Pt{L(K) < 9" < U(K) | 9} > min (?t{L(K) < 9' < U(K) \ 9'}, Pr{L(K) < 9" < U{K) \ 9"}) 

for any 9 € (9', 9"), where the lower bound in the right side is no greater than mm(Pr{L(K) < 
9' < U(K) | 9'}, Pi{L{K) < 9" < U{K) \ 9"}). Hence, 

inf Pr{i(A") < 9 < U{K) \ 6} > min (Pi{L(K) < 6' < U{K) \ 9'}, Pr{L{K) < 9" < U(K) \ 9"}) . 

ee[e',e"] 

Now we show that the equality must hold. For simplicity of notations, let the left and right sides 
in the above inequality be denoted by 7 and p respectively. Suppose 7 > p. Then, 

Px{L(K) <9< U{K) I 9} > W € [9',9"], 

which implies that 

Pt{L{K) < 9' < U(K) I 9} = Pv{L(K) < 9" < U(K) \ 9} > G (9', 9"). 

Recalling that both Px{L{K) < 9' < U{K) \ 9} and Pr{L(K) < 9" < U{K) \ 9} are continuous 
for 9 G [a, b], we have 

p = min (Pt{L(K) < 9' < U(K) | 9'}, Pt{L(K) < 9" < U{K) j 9"}) > 

leading p > 7, which contradicts 7 > p. Therefore, it must be true that 7 = p. That is, 

inf Pt{L{K) <6 < U{K) I 6} = min (Pr{i(A") < 6' < U{K) \ 9'}, Pr{L{K) < 9" < U(K) \ 6"}) . 

It follows that the infimum of Pr{L(K) < 9 < U{K) \ 9} with respect to 9 G [a, b] equals the 
minimum of the set {C(a), C{b)} U {C v {6) : 9 G £lu U £ L } U {C L {9) : 9 G J3l U Let 

£' = £u n g L , % = £u\ £' L = £ L \ 

Then, 

{C v {9) : 9 G J2jj U J2l} U {C l {9) :9e£uU£ L } 
= {C v {9) :9£%}U {Cu{9) : 9 G U {C L {9) : 9 G B' L ) U {C L {9) : 9 G JS'} 

u{Cu(9):9eJ2' L }u{C L (9):9e%}. (7) 
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For 9 G we have < Pv{L(K) = 9 < U{K) \ 9} < Pr{L{K) = 9 \ 9} = and 

Cu(9) - C L {9) = Pr{L(K) <9< U(K) \ 9} - Pr{L(K) <9< U(K) j 9} 
= Pv{L(K) = 9 < U(K) j 9} - Pr{L(K) <9 = U(K) \ 9} 
= - Pt{L(K) <9 = U{K) | 9} < 0, 

which implies that 

min {Cjj{9) : 9 G %} < min {C L {9) : 9 G %} . (8) 
For 9 G £' L , we have < Pv{L(K) <9 = U(K) \ 9} < Pr{U(K) = 9 \ 9} = and 

Cu(0) - C L {9) = Pr{L(K) = 9< U{K) j 9} - Pr{L(K) <9 = U(K) | 9} 
= Pi{L{K) = 9 < U(K) \9}>0, 

which implies that 

min {C v (9) : 9 G & L } > min {C L {9) : 9 G £' L } . (9) 
Combing ©, (HD and ([9]), we have 

min {Cu(9) : 9 G £?u U £ L } U {C L {9) : 9 G JSu U B L } 
= min {Cjj{9) : 9 G %} U (0/(0) : 9 G =2'} U (C L (#) : G ^} U {C L {9) : 9 G £'} 
= min {Ccr(0) : 9 G U {C L (fl) : 9 G , 

which implies that the minimum of the set {C(a), C(b)} U {Cu{9) : 9 G B v U U {C L (6>) : (9 G 
^ U B L } equals the minimum of {C(a), C(6)} U {0/(0) : 9 G U {C L (0) : G ^l}- This 
proves statement (II) of Theorem 7. 

Clearly, statement (III) of Theorem 7 is already justified in the course of proving statements 
(I) and (II). This concludes the proof of Theorem 7. 



B Proof of Theorem 8 

For the simplicity of notations, define 

,, m! v . if0<z<m, 

z\(m—z)l — — ' 

if z < or z > m 

for non-negative integer m and arbitrary integer z. We now establish some preliminary results. 

Lemma 2 Let < M < N. Define T(k, M, N, n) = ( J ]f ) (^iJJ^ 1 ) / Q- Then, Pv{K < k 
M} — Pr{K < k | M + 1} = T(k, M, N, n) for any integer k. 
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Proof. We first show the equation for < k < M. We perform induction on k. For k = 0, we 
have 

Pv{K <k\M}- Pr{K < k \ M + 1} = Pv{K = | M} - Pv{K = | M + 1} 



(D G) 

/iV-M-l\ 
I n-1 / 

o 

/M\ /JV-Af-l\ 
I A n-0-1 / 



where (|10p follows from the fact that, for non-negative integer m, 



T(0, Af, AT, n) 



z + \z ) \z+l 

for any integer z. 

Now suppose the lemma is true for k — 1 with 1 < A; < M, i.e. 



Pr{ET < jfe - 1 | M} - Pr{K <k-l\M + l} = — 



Then, 

Pt{K <k\M}- Pr{K < k | M + 1} = Pr{it < fc - 1 | M} - Pr{K < k - 1 | M + 1} 

(M+l\ (N-M-l\ 
\ k I \ ri-fc / 



o 


o 






/ M \ (N-M-l\ 
U-lA n-k I 


(M\ (N-AI\ 
\k)\ n-k I 


o 


o 






\k)\ n-k I 


\ k )\ n-k ) 


o 


[ o 






(M\ (N-AI\ 
\k)\ n-k ) 


V A. / \ it — fin / 


o 


o 






\ k) \ n-k-1 I 




o 









a 

M \ (N-M-l\ 



(10) 



l) \ n-k ) 



o 

(12) 
(13) 

where (|12p and (|13p follows from (jlip . Therefore, we have shown the lemma for < k < M. 

For k > M, we have Pr{K < k \ M} = Pr{K <fc|M + l} = l and T(k, M, N, n) = 0. For 
k < 0, we have Pr{i<C < k \ M} = Pi{K < k \ M + 1} = and T(k, M, N, n) = 0. Thus, the 
lemma is true for any integer k. 

□ 

Lemma 3 Let 1 < M < N and k < I. Then, 

Pv{k < K < I | M} - Pr{k < K < I | M - 1} = T(k - 1, M - 1, N, n) - T(Z, M - 1, N, n). 
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Proof. To show the lemma, it suffices to consider 6 cases as follows. 

Case (i): < n < k < I. In this case, Pv{k < K < I \ M} = Pi{k <K<I|M-1} = and 
T(k — l,M — l,N,n) = T(l, M - 1, N, n) = 0. 

Case (ii): k < I < < n. In this case, Pi{k < K < I \ M} = Pr{k < K < I \ M - 1} = and 
T(k — 1,M — 1,N, n) = T(l, M - 1, N, n) = 0. 

Case (iii): k < < n < I. In this case, Pr{k <K <l\M} = Pi{k < K < I | M - 1} = 1 and 
T(k — l,M — l,N,n) = T(l, M - 1, N, n) = 0. 

Case (iv): k < < I < n. In this case, T{k — 1, M — 1, N, n) = and, by Lemma [21 

Pr{fc < K < I | M} - Pr{fc < K < I \ M - 1} = Pr{A" < I \ M} - Pr{A' < I \ M - 1} 

= T(k- 1,M- l,N,n) -T(l,M-l,N,n). 

Case (v): < k < n < I. In this case, T(l, M — 1, N, n) =0 and, by LemmaEl 

Pr{fc < K < I | M} - Pr{fc < K < I \ M - 1} = Pr{A" < k \ M - 1} - Pr{K < k \ M} 

= T(k- 1,M- l,N,n) -T(l,M-l,N,n). 

Case (vi): < k < I < n. In this case, by Lemma [2j 

Prjfc < K < I | M} - Pr{k < K < I \ M - 1} 
= [Pr{A" <l\M}- Px{K < k \ M}\ - [Pr{K < I \ M - 1} - Px{K <k\M — 1}] 
= [Pr{A <l\M}- Px{K <l\M-l}] - [Pr{K < k \ M} - Pr{K < k \ M - 1}] 
= T(k- l,M-l,N,n) -T(l,M- l,N,n). 



□ 



Lemma 4 Let I > and k < n. Then, 
M<1+ 

— n— 1 



nM 



> I for M > 1 + 



71-1 



and 



nM 
N+l 



< k — 1 /or 



Proof. To show the first part of the lemma, observe that (N + 1 — n)l > 0, by which we can 



show %1±>(N+ 1)1. Hence, n [ 1 + 



n-l 



> ^{ > (iV + 1)L That is, 



JV+l 



1 + 



Nl 
n-l 



> I. 



It follows that 

nM 



N+l 



N+l 

> I for M > 1 + 



1 + 



n — 1 



> L Since the floor function is non-decreasing, we have 



Nl 
n-l 



To prove the second part of the lemma, note that (N + 1 — n){n — k) > 0, from which we 

< k, 



can deduce 1 + < ^"t • Hence, 1 + 



AT(fc-l) 
n-l 



(JV+l)fc . 
^ n ' i-e! -> AT+1 



1 + 



N{k-1) 
n-l 



leading to 



Af+l 



1 + 



N(k-l) 
n-l 



< k — 1. Since the floor function is non-decreasing, we have 



nM 
N+l 



<fc-lforM<l + 



jV(fc-l) 
n-l 



□ 
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Lemma 5 Let < r < re. Then, the following statements hold true. 
(I) 

T(r — 1, M - 1, N,n) < T(r, M — 1,N, re) for 1 < r < 



(II) 



T(r + 1,M- l,N,n) < T(r,M- l,N,n) for 



reM 
N + l 



T(r, M - 2, N, n) < T(r, M — 1,N, n) for 1<M<1 

Nr 

T(r, M, N, n) < T(r, M — 1,N, n) for 1 + 



nM 
N + l_ 

< r < n — 1. 
Nr 



n — 1 



re — 1 
< M < N. 



Proof. To show statement (I), note that T(r, M — 1, N,n) = for min(M — 1, re — 1) < r < re. 
Our calculation shows that 

T(r - 1,M - 1,N, re) r iV — M + 1 — (re - r) „ reM 

< 1 for 1 < r < 



and 



T(r,M - 1,N, re) M-r re-r 



T(r - 1,M - 1,N, re) reM 

— H ' ! — t-^ > 1 for < r < min(M - 1, re - 1) 

T(r,M - l,N,n) Ar 1 1 ~ v 



iV+1 



iV + 1 

To show statement (II), note that T(r,M - 1, iV, re) = for 1 < M < r + 1, and T(r, M 
1, iV, re) > T(r, M — 2, N,n) = for M = r + 1. Direct computation shows that 

T(r, M — 1, N,n) M - 1 N - M + 2 - (n - r) r Nr 

■ > 1 for r + 1 < M < 1 1 



T{r, M - 2, N, n) M — 1 — r N-M+l 



n — 1 ' 



and 



T(r,M — l.N.n) „ Nr 

' < 1 for 1 + r <M <N. 



T(r, M-2,N, re) 



re — 1 



□ 



Lemma 6 Let < £ <U < N . Then, for any integers k and I, Pr{/c < K < I \ M} is unimodal 
with respect to M for C < M <U. 

Proof. Clearly, the lemma is trivially true if k > I. Hence, to show the lemma, it suffices to 
consider 6 cases as follows. 

Case (i): < re < k < I. In this case, Pr{fc <K<1\ M} = for any M € [C,U]. 

Case (ii): k < I < < re. In this case, Pr{k < K < I | M} = for any M £ [£,U\. 

Case (hi): k < < re < I. In this case, Pv{k < K <l \ M} = 1 for any M £ [C,U]. 

Case (iv): k < < I < re. In this case, Pr{fc < K < I \ M} = Pr{K < I \ M} is non- increasing 
with respect to M € [C,K\ as can be seen from Lemma [21 

Case (v): < k < n < I. In this case, Pr{/c < K < I \ M} = 1 - Pi{K < k | M} is 
non-decreasing with respect to M € [£,14] as can be seen from Lemma [2j 
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Clearly, the lemma is true for the above five cases. 

Case (vi): < k < I < n. Define A{k,l,M,N,n) = Pr{k < K < I \ M} - Pr{fc < K < I \ M- 1}. 
By Lemmal A(k, I, M, N, n) = T(k — 1,M — 1, N, n) - T(l, M - 1, N, n). 



Invoking Lemma U for M > 1 + 



Nl 

n-1 



, we have that 



nM 
N+l 



> I and thus, by statement (I) of 



Lemma El T(r, M — 1, N, n) is non-decreasing with respect to r < I. Consequently, T{k — 1, M 
1, N, n) < T(l, M - 1, N, n), leading to A(k, I, M, N, n) < for M > 1 + 



Similarly, applying Lemma HI for M < 1 + 



N{k-l) 
n-1 



we have that 



Nl 
n-1 

nM 
N+l 



< k — 1 and thus, 



by statement (I) of Lemma [5l T(r,M — l,N,n) is non-increasing with respect to r > k — 1. 
Consequently, T(k - 1,M - l,N,n) > T(/,M - l,N,n), leading to A(k,l,M,N,n) > for 
M < 1 + . 

Nl 



By statement (II) of Lemma [SJ for 1 + 



AT(fc-l) 
n-1 



< M < 1 + 



n-1 

1,M - 



we have that T(l, M — 
l,N,n) is non- increasing 



l,N,n) is non-decreasing with respect to M and that T(k 
with respect to M. It follows that A(k,l,M,N,n) is non-increasing with respect to M in this 



JV(fc-l) 
n-1 



< M* < 1 + 



jVj 
n-1 



and that 



range. Therefore, there exists an integer M* such that 1 + 
A{k,l,M,N,n) > for < M < M* , and A(k,l,M,N,n) < for M* < M < N. This implies 
that Pr{/c < K < I \ M} is non-decreasing for < M < M* and non-increasing for M* < M < N . 
This concludes the proof of the lemma. 

□ 



Lemma 7 Let < M < N. Then, Pr{g <K<h + l\ M + l}> Pr{g < K < h \ M} for any 
integers g and h. 

Proof. Clearly, the lemma is trivially true if g > h. Hence, to show the lemma, it suffices to 
consider the case g < h. Note that, by Lemma [31 



Pi{g < K < h + 1 | M + 1} - Vi{g <K <h\M} 
M + l\(N-M-l\ UN 
h+1 )\ n- h-1 ) I \n 
M + 1\(N- M — 1\ UN 
h + 1 J\n-h-l ) / \n 
M + 1\(N-M-1\ (M\ fN — M — l 
h + l)\n-h-l ) ~ \h J\ n-h-1 
M \ fN — M — 1\ I (N 



+ Pv{g <K<h\M + l}- Pr{g < K < h \ M} 
+ T(g - 1, M, N, n) - T{h, M, N, n) 



N 



n 



+ T(g-l,M,N,n) 



K h + 1/ \ n — h — 1 j 
where the last equality follows from (jlll) . 



+ T(g-l,M,N,n)>0, 



□ 



Lemma 8 Let < M < N. Then, Pi{g -l<K<h\M-l}> Pr{g < K < h \ M} for any 
integers g and h. 
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Proof. Clearly, the lemma is trivially true if g > h. Hence, to show the lemma, it suffices to 
consider the case g < h. Note that, by Lemma [3l 



Pr{# -l<K<h\M-l}- Yi{g < K < h | M} 



M - 1 

9-1 
M - 1 

9-1 
M - 1 

9-1 
M - 1 

5-1 



N-M + l 
n — g + 1 
N-M + l 
n — g + 1 
N-M + l 
n- g + 1 
N-M 
n — g + 1 



N 

N 

n 
M 
9 - 



+ Pv{g < K < h | M - 1} - Pr{# < K < h \ M} 
+ T(h, M - 1, N, n) - T(g — 1,M — 1,N, n) 



1 



N-M 
n-g 



+ T(h,M - l,N,n) 



+ T(h,M- l,N,n) > 0, 



where the last equality follows from ([TT 



□ 



Lemma 9 Suppose that {M 1 < L{K) < M"} = {M 1 < U(K) < M"} = 0. Then, Pi{L(K) < 
M < U(K) | M} is unimodal with respect to M for M' < M < M" . 

Proof. First, we shall show the following facts: 

(i) If {L{K) = M'} = 0, then {L{K) < M} = {L{K) < M'} = {L(K) < M"} for M' <M< 
M". 

(ii) If {L(K) = M'} ^ 0, then {L(K) < M} = {L(K) < M'} = {L(K) < M"} for M' < M < 
M" . 

(iii) If {U(K) = M"} = 0, then {U(K) > M} = {U{K) > A/'} = {U(K) > M"} for 
M' <M< M" . 

(iv) If {U{K) = M"} + 0, then {U(K) > M} = {U{K) > M'} = {U{K) > M"} for 
M' < M < M". 

To show statement (i), making use of {L(K) = M'} = {M' < L(K) < M"} = 0, we have 
{M' < L(K) < M} = {M' < L(K) < M} C {M' < L(K) < M"} = and {L(K) < M} = 
{L{K) < M'} U {M' < L(K) < M} = {L(K) < M'} for M' < M < M" . On the other hand, 
{L(K) < M} = {L(K) < M"} \ {M < L{K) < M"} = {L{K) < M"} for M' < M < M" . 

To show statement (ii), making use of {M 1 < L(K) < M"} = 0, we have {M 1 < L(K) < 
M} C {M' < L(K) < M"} = and {L(K) < M} = {L(K) < M'} U {M 1 < L{K) < M} = 
{L(K) < M'} for M' < M < M" . On the other hand, {L(K) < M} = {L(K) < M"} \ {M < 
L(K) < M"} = {L(K) < M"} for M' < M < M" . 

To show statement (iii), using {U(K) = M"} = {M' < U{K) < M"} = 0, we have {M' < 
U(K) < M} C {A/' < U(K) < M"} = and {U(K) > M} = {U(K) > M'} \ {AT' < U(K) < 
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M} = {U(K) > M'} for M' < M < M". On the other hand, {U(K) > M} = {U{K) > 
M"} U {M < U(K) < M"} = {U(K) > A I"} for M' < M < M" . 

To show statement (iv), note that {U(K) > M} = {U(K) > M'} for M' < M < M" . On 
the other hand, {U(K) > M} = {U(K) > M"} U {M < U(K) < M"} = {U(K) > M"} for 
M' < M < M". 

Now, to show the lemma, it suffices to consider four cases as follows. 

Case (i): {L(K) = M'} = 0, {U{K) = M"} = 0. 
Case (ii): {L(K) = M'} = 0, {U{K) = M"} ^ 0. 
Case (iii): {L(K) = M'} f 0, {U(K) = M"} = 0. 
Case (iv): {L(K) = M'} + 0, {U(K) = M"} + 0. 

In Case (i), making use of facts (i) and (iii), we have {L{K) < M < U(K)} = {L(K) < 
M' < U(K)} = {L(K) < M" < U(K)} for M' < M < M" . Invoking Lemma El we have that 
Pr{L(K) < M < U(K) \ M} is unimodal with respect to M for M' <M< M" . 

In Case (ii), making use of facts (i) and (iv), we have {L(K) < M < U(K)} = {L(K) < 
M' < U{K)} = {L(K) < M" < U(K)} for M' < M < M" . Invoking Lemma we have 
that Pi{L(K) < M < U(K) \ M} is unimodal with respect to M for M' < M < M" . Since 
{M" = U(K)} / and U(K) is monotonically increasing, we have {M" < U(K)} = {K > k} 
and {M" < U(K)} = {K > k + 1}, where k = mm{k : U{k) > M"} <k = max{k : U(k) < M"}. 
Therefore, as a result of Lemma [HJ 

Pr{L(K) < M" < U(K) \ M" - 1} > Pr{L(K) < M" < U(K) | M"}. (14) 

It follows that Pr{L(K) < M < U(K) \ M} is unimodal with respect to M for M' < M < M" . 

In Case (iii), making use of facts (ii) and (iii), we have {L(K) < M < U(K)} = {L(K) < 
M' < U(K)} = {L(K) < M" < U(K)} for M' < M < M". Invoking Lemma El we have 
that Pr{L(K) < M < U(K) \ M} is unimodal with respect to M for M' < M < M" . Since 
{M' = L(K)} 7^ and L(K) is monotonically increasing, we have {M' > L(K)} = {K < k} and 
{M' > L(K)} = {K < k- 1}, where k = min{k : L(k) > M'} <k = max{k : L(k) < M'}. 
Therefore, as a result of Lemma [HJ 

Pr{L(K) < M' < U(K) \ M'} < Pv{L(K) < M' < U(K) | M' + 1}. (15) 

It follows that Pt{L(K) < M < U(K) \ M} is unimodal with respect to M for M' <M< M" . 

In Case (iv), making use of facts (ii) and (iv), we have {L(K) < M < U(K)} = {L(K) < 
M' < U(K)} = {L(K) < M" < U(K)} for M' < M < M" . Invoking Lemma El we have that 
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Pr{L(K) < M < U(K) \ M} is unimodal with respect to M for M' < M < M" . Recalling 
(fT4|) and (US]), we have that Vi{L{K) < M < U(K) \ M} is unimodal with respect to M for 
M' <M< M" . 

□ 

Finally, we are in a position to prove the theorem. Let M' < M" be two consecutive distinct 
elements of I UL . Then, {W < L(K) < M"} = {W < U(K) < M"} = 0. By LemmaEl we have 
that Pr{L(K) < M < U(K) \ M} is unimodal with respect to M for M' < M < M" . Since this 
argument holds for any consecutive distinct elements of the set Iul, Theorem 8 is established. 
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